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ABSTRACT 

We calibrate the galaxy cluster mass - temperature relation using the temperature profiles of intr- 
acluster gas observed with ASCA (for hot clusters) and ROSAT (for cool groups). Our sample con- 
sists of apparently relaxed clusters for which the total masses are derived assuming hydrostatic equi- 
librium. The sample provides data on cluster X-ray emission-weighted cooling flow-corrected temper- 
atures and total masses up to riooo- The resulting M-T scaling in the 1-10 keV temperature range is 

Miooo = 1-23 ± 0.20 lO^s h^^ Mq ( j with 90% confidence errors, or significantly (99.99% 

confidence) steeper than the self-similar relation M oc T^/^. For any given temperature, our measured 
mass values are significantly smaller compared to the simulation results of Evrard et al. (1996) that 
are frequently used for mass-temperature scaling. The higher-temperature subsample (kT > 4 keV) 
is consistent with M cx T^/^, allowing the possibility that the self-similar scaling breaks down at low 
temperatures, perhaps due to heating by supernovae that is more important for low-temperature groups 
and galaxies as suggested by earlier works. 

Subject headings: cosmology: observations ~ dark matter - galaxies: clusters: intergalactic medium - 
X-rays: galaxies 



1. INTRODUCTION 

Galaxy clusters are the largest gravitationally bound 
objects in the universe, and thereby provide information 
on cosmic structure formation. The mass distribution of 
virialized objects can be predicted for different cosmolo- 
gies and different initial density fluctuation spectra. By 
comparing such predictions to the observed cluster mass 
function, one can constrain cosmological parameters. Cos- 
mological parameters most strongly affect mass function 
predictions for large masses that correspond to the cluster 
scale, and therefore clusters are of special importance. 

Accurate measurements of the cluster total mass, domi- 
nated by dark matter, are challenging at present and pos- 
sible only for a limited number of clusters. For this reason 
there is currently insufficient data for a direct derivation 
of the mass function. A more practical way of determining 
the mass function is to observe the distribution of readily 
available average cluster gas temperatures and to convert 
this to a mass function, taking advantage of the tight mass 
- temperature correlation predicted by hydrodynamic clus- 
ter formation simulations (e.g. Evrard, Metzler & Navarro 
1996). The mass and temperature are predicted to scale 
as M oc Although different simulations and ob- 

servations are in general qualitative agreement, there are 
significant disagreements on the details of the gas temper- 
ature profiles (Frenk et al. 1999, Markevitch et al. 1998). 
Therefore, this relation needs observational confirmation 
and calibration. 

In addition to being useful for providing a link between 
cosmological predictions and observations, the M-T rela- 
tion is also interesting in itself, because any deviations 
from the predicted self-similar scaling of Al oc T^/^ would 
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indicate that additional physical processes are at play than 
gravity alone. The detailed behaviour of the M-T rela- 
tion provides information about the process of cluster for- 
mation and energy input into the gas. One particularly 
important energy input source is preheating of the inter- 
galactic gas by early supernova driven galactic winds (e.g. 
David, Forman & Jones (1991), Evrard & Henry (1991), 
Kaiser (1991), Loewenstein & Muskotzky (1996)). 

Steps toward calibrating the M-T relation with obser- 
vations have been made by measuring cluster masses us- 
ing gravitational lensing (Hjorth, Oukbir & van Kampen 
1998) and the hydr ostatic equilibrium approa ch assum- 
ing isothermal gas (Neumann & Arnaud 1999). Horner, 
Mushotzky & Scharf (1999) derived M-T relation using 
several independent methods: virial theorem for cluster 
galaxies and hydrostatic equilibrium applied to the X- 
ray emitting gas, both assuming isothermality and using 
published mass values derived from measured temperature 
profiles. 

As in the latter works, in this paper, we derive a clus- 
ter mass - temperature relation under the assumption of 
hydrostatic equilibrium. We use the published total mass 
profiles and X-ray emission- weighted temperatures derived 
from temperature profiles of hot clusters measured with 
ASCA. For cooler groups, using the published ROSAT 
cluster temperature profiles, we compute the correspond- 
ing temperature values, and the total mass profiles in cases 
where masses are not published in sufficient detail. We in- 
corporate several improvements over the work of Horner 
et al. (1999). Our ASCA cluster sample is homogeneous, 
the temperature profiles are all determined using the same 
method that accounts for the ASCA PSF (Markevitch et 
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al. 1998) so that the resulting mass values and their er- 
rors are directly comparable. The Horner et al. (1999) 
sample contains mass profiles derived both with ASCA 
PSF correction (e.g. A2256 and A2029) as weU as with- 
out it (e.g. A496 and A2199). For A496 and A2199 the 
PSF correction is not large for the central pointing used 
in Horner et al. (1999), but we measure the temperature 
profile to a larger radius using offset pointings where the 
PSF correction is significant. Another important differ- 
ence is that Horner et al. (1999) extrapolate mass pro- 
files to an overdensity radius r2oo assuming pdark oc r~^-^, 
whereas we use the measured mass profiles up to riooo ■ At 
this radius relaxed clusters are unlikely to experience sig- 
nificant residual turbulence that would violate the hydro- 
static equilibrium assumption (e.g. Evrard et al. 1996). 
Our sample contains only clusters without signs of dis- 
turbance. Thirdly, we combine our ASCA sample with a 
low-temperature subsample of galaxy groups and galax- 
ies, whose temperature profiles have been measured with 
ROSAT, in order to study the M-T relation over a tem- 
perature range of 1 - 10 keV. 



We use Ho = 50/i5okms ^ Mpc" 
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1 and report 



90% confidence intervals throughout the paper, except 
where stated otherwise. 

2. DATA 
2.1. The sample 

The sample consists of 9 clusters, groups and galaxies 
with pubhshed (see Section 2.3) relatively accurate spa- 
tially resolved temperature profiles measured with ASCA 
(the 6 hotter ones) and ROSAT (the 3 cooler ones) up 
to radii of overdensity of 1000 - 500, comparable to the 
virial radii. Hereafter we use the usual notation and 
Mi where is a radius of overdensity (the mean interior 
density with respect to the critical density) of i, and Mi is 
the mass within that radius. For our sample, we required 
the objects to be apparently relaxed, with no substructure 
or deviation from azimuthal symmetry in ROSAT PSPC 
images, or in the ASCA temperature maps in Markevitch 
et al. (1998). Hence, the effects of bulk motions, and 
the deviation from hydrostatic equilibrium are minimized. 
The requirement of apparent relaxation limits severely the 
number of suitable cluster candidates for our analysis. 

Most relaxed clusters have strong cooling flows which 
complicate the ASCA spatially resolved analysis. Since 
the ASCA PSF has a half-power diameter comparable to 
the angular size of a typical cooling flow in a nearby clus- 
ter, the cooling flow parameters can not be constrained ad- 
equately. However, when modeling the temperature struc- 
ture of the cluster, the (uncertain) cooling flow model must 
be included as a component. The wide energy-dependent 
PSF scattering of the photons from a strong cooling flow 
therefore increase the uncertainty of the non-cooling-flow 
model temperatures even at a large radius. Therefore we 
mostly selected relaxed clusters with only moderate cool- 
ing flows (and therefore accurate temperature proflles), 
which further limited the number of suitable sample clus- 
ters. 

We seek relaxed systems which span the 1-10 keV tem- 
perature range. The clusters with the lowest temperatures 
are also the faintest and accurate spectroscopy is possible 
only for the nearest ones. The obvious problem is that 



in nearby clusters the virial radius corresponds to a large 
angular distance, usually beyond the ASCA and ROSAT 
field of view. This furthermore limits the number of avail- 
able objects at low temperatures (kT < 4 kcV). There 
are several nearby cool groups and galaxies, for which the 
analysis with the needed accuracy is feasible with ROSAT, 
with a slight extrapolation of the mass profiles (see 2.4.). 

Our search in the archives and literature resulted in 
the following sample of objects suitable for our analysis: 
A496, A2199, A401, A3571, A2256 and A2029 studied 
with ASCA, and NGC5044, NGC507 and HCG62 with 
ROSAT. 



2.2. Th 



e mass errors 



Total masses are determined assuming hydrostatic equi- 
librium using the measured gas temperature and density 
profiles. For some clusters the formal uncertainty of the 
resulting mass values is small. However, hydrodynamic 
simulations show that the systematic uncertainties inher- 
ent in the hydrostatic mass determination method, such as 
deviations from spherical symmetry or hydrostatic equi- 
librium due to incomplete thermalization of the gas, will 
lead to about a 10 - 30% uncertainty in the calculated 
total masses (e.g. Evrard et al. 1996; Schindler 1996; 
Roettiger, Burns & Loken 1996). Since these simulations 
include also mergers, whereas our sample contains only 
relaxed clusters, our uncertainties would be towards the 
low end of the above interval. However, to be conserva- 
tive, if any formal mass error from the literature is smaller 
than 20% of the mass value we use a 20% uncertainty in- 
stead. The errors on the emission- weighted temperatures 
are much smaller than the mass errors. For example, at 
the most interesting radius riooo (see below) most of the 
relative temperature errors are smaller than 0.2 times the 
relative Afiooo errors wherefore we ignore the temperature 
errors while deriving the M-T relation. 

2.3. The ASCA subsample 

For the hotter (> 4 keV) clusters (A496, A2199, A401, 
A3571, A2256, A2029) we use the emission weighted 
cooling-flow corrected ASCA temperatures from Marke- 
vitch et al. (1998). For these clusters we use the published 
total mass profiles obtained using observed ASCA temper- 
ature proflles and ROSAT s urface brightness profiles . The 
masses of A496 and A2199 ( [Markevitch et al. 19991 ) have 
been derived by modeling the temperature profile with a 
polytropic form an d the mass for A2029 (Sarazin, Wise 
fc Markevitch 1998| ) has been derived by modeling a tem- 
perature profile as a linear function of the radius. For 
A401 (Nevalainen, Markevitch & Forman 1999a), A3571 
(Nevalainen, Markevitch & Forman 1999b) and A2256 
(Markevitch & Vikhlinin 1997) the masses have been ob- 
tained modelling the dark matter density with different 
functional forms (see above papers for the details), com- 
puting the corresponding temperature profiles and fitting 
the dark matter distribution parameters. For all these 
clusters the temperature profile is determined relatively 
accurately out to r^oQ. All these clusters exhibit a tem- 
perature decline with radius, which results in smaller total 
mass values at our radii of interest, compared to the usual 
isothermal analysis. The "universal" NFW dark matter 
profile (Navarro, Frenk & White 1997) provides a good 
description of these profiles over a large range of radii. 
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2.4. The ROSAT subsample 

For the cooler objects, we used the pubhshed ROSAT 
temperature and surface brightness profiles to compute the 
emission weighed mean temperatures outside the cooling 
flow regions, to be directly comparable to the cooling flow- 
corrected ASCA temperatures for hotter clusters. We wish 
to do all the analysis at rsoo, but the ROSAT PSPC data 
do not extend to such large radii. We extrapolate the mass 
profiles beyond the radii of measured temperatures (except 
in the case of NGC507, where the temperature profile ex- 
tends sufficiently far), but not beyond the radii to which 
the cluster surface brightness is significantly detected. For 
the group HCG62 and the galaxy NGC507 this allows us 
to reach riooo, but for the NGC5044 group we reach only 
''1500 (see Figure 1). While for the hotter clusters we use 
the published mass and temperature values, ROSAT data 
need some additional analysis, as described in detail below. 

2.4.1. Group HCG62 

Using the ROSAT PSPC temperature profile from Pon- 
man & Bertram (1993) we obtain the X-ray emission- 
weighted temperature (T) = 1.16 ± 0.08 keV beyond the 
cooling radius of r = 2.4'. For comparison, a published 
single component temperature of ASCA data f or HCG62 
gives a smaller value of (T) = 1.05 ± 0.02 keV ( [Fukazawcj 
et al. 1998|). However, in the ASCA analysis the central 



5Q Mpc was excluded, which corresponds to 6.4' for 



0.15/1, 

HCG62. Thus the ASCA value integrates over the outer 
parts of HCG62. In that region, the dechning ROSAT 
temperature profile is consistent with the above ASCA 
value. Therefore, we use the value derived from the profile 
of Ponman & Bertram (1993). 

To derive the mass profile, we fitted the above ROSAT 
temperature profile with a poly tropic form (T(r) oc 
Pgas{f')'^~^), fixing the gas density using the results of the 
ROSAT surface brightness analysis (Ponman & Bertram 
1993), and found the best fit with 7 = 1.09t^;i|. For this 
fit, we excluded the cooling fiow region r < 2.4'. The total 
mass for a polytropic temperature profile inside a radius 



is given by 



Mtotir) = 3.70 X 10^3 jvf T(r) a, 

1 -I- a;- 



a;2 0.60' 



(1) 



where /i is the mean molecular weight, /? and are the 
slope parameter and core radius of the gas density profile, 
T is expressed in keV and ax in Mpc. We propagated the 
errors of 7 and T to the total mass. 

2.4.2. Galaxy NGC507 

Kim & Fabbiano (1995) obtained a temperature profile 
for this galaxy with ROSAT PSPC. Excluding the data 
from the cooling flow area (r < 2.0'), we obtain a mean 
temperature of (T) = 1.09 ± 0.03 keV. The corresponding 
value obtain ed by ASCA is significa ntly higher, {T)— 1.26 
± 0.07 keV (Fukazawa et al. 1998). We are not strongly 
concerned about this disagreement, because the ROSAT 
and ASCA energy bands are very different and there may 
be differences in the average temperature due to, for ex- 
ample, nonisothermality and a possible nonthermal com- 
ponent. We assume that for low temperature systems ( 
(T) ~ 1 keV), ROSAT gives a more accurate value, and 



use the above ROSAT value in our analysis. We show 
below that if we use the higher ASCA temperature value 
instead, our conclusions will only be strengthened. 

The total mass profile is directly taken from Kim & Fab- 
biano (1995), who used a polyt ropic model for the temper- 
ature profile, as described in (2.4.1) above. The errors of 
temperature and the polytropic index 7 are incorporated 
in the mass values. 

2.4.3. Group NGG50U 

Using the ROSAT temperature profile of David et al. 
(1994), we obtained the emission weighed temperature {T) 
= 1.09 ± 0.03 keV, beyond the cooling radius r = 4.0'. 
ROSAT temperatures beyond O.lbh^Q Mpc are consistent 
with the corresponding A SCA value ((T) = 1.07 ± 0.01 
keV [Fukazawa et al. 19981 ). 

We use the analytical form for the best fit total mass 
given in David et al. (1994), obtained using a power-law 
model for the temperature profile. We propagate the er- 
rors of temperature and the exponent of the temperature 
profile to the total mass errors. Because the X-ray emis- 
sion is not detected beyond risoo for this nearest group, we 
will not use this group in our analysis beyond that radius. 

2.4.4. NFW mass profile 

The NFW model for the dark matter density profile de- 
scribes well the mass distribution in hot ASCA clusters in 
our sample (s ee e.g. Markevitch fc Vikh linin 1997, Sarazin 



et al. 1998, Markevitch et al. 1999, Nevalainen et al 



1999a, Nevalainen et al. 1999b). I'urthermore, a similar 
result was found for Coma in an optical study of galaxies 
(Geller, Diaferio & Kurtz 1999). Therefore, we compared 
the NFW model with the derived mass profiles for the cool 
groups and galaxies HCG62, NGC507 and NGC5044. We 
fixed the gas density profile to ROSAT values and modeled 
the dark matter density distribution with an NFW profile 
selecting its parameters so that gas plus dark matter ap- 
proximate the derived total mass profiles. For all three 
systems, the NFW profile provides a good description of 
the dark matter distribution over a range of interesting 
radii (see Figure 1). Such an underlying universal dark 
matter profile is con sistent with the observed s imilarity of 
gas density profiles (Neumann & Arnaud 199£, Vikhlinin, 
Forman & Jones 1999) and temperature profiles (Marke- 
vitch et al. 1998) in different clusters, when scaled to 
physical radii. 

3. RESULTS 

For each object, using its measured total mass profile, 
we computed the overdensity, or the mean interior density 
in units of the critical density as a function of the radius: 
Mtotir) /i^irr^Pc), where pc = 3H^{1 + z)^ /8ttG. We then 
calculated the masses within several radii of fixed over- 
density (see Table 1). To remove the different amounts 
of evolution in our sample temperatures due to different 
values of 2, we scaled the observed mean temperatures to 
z = according to the predicted scaling Tgas{z) oc (1 -I- z) 
for a given mass (e.g. Eke, Cole & Frenk 1996). To be 
consistent with this, we should also compute the overden- 
sities at z = 0, but this would require a further assumption 
that we observe the clusters just after their collapse. How- 
ever, due to low z values in our sample, our results do not 
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change significantly, whether wc compute the overdensities 
at z = or at the observed rcdshifts, and we report our 
results using the latter method. 

Figure 2 shows that Miooo is strongly correlated with 
(T), the X-ray emission weighted cooling- flow corrected 
temperature. We fit the masses with a function 

where i denotes different values of overdensity. 

We are able to obtain interesting constraints for the 
above relation for radii up to riooo- We exclude the low 
temperature data from the fits at larger radii, because the 
gas in the nearby ROSAT objects is not detected to such 
low overdensities. At riooo we get an acceptable fit with 
mooo = l-23l°:?^ X W Mq, aiooo = l-79t°;}t 
= 6.4 for 6 d.o.f. The self-similar prediction a = 3/2 
differs by 3.7 a. If we use the higher ASCA tempera- 
ture for NGC507 (see 2.4.2) the relation steepens by a few 
per cent, and the difference with the self-similar slope in- 
creases slightly. If we were to include the data of NGC5044 
in the fit at riooo (see 2.4.3.), the results would not change 
significantly since that data point is consistent with the 
ones of HCG62 and NGC507, but the constraint would 
improve. Therefore, our choices of using the lower temper- 
ature for NGC507 and excluding NGC5044 data at riooo 
in all our fits are conservative. Fixing the slope aiooo to 
1.5 the best fit gives = 29.4 for 7 d.o.f., thus aiooo = 
3/2 is ruled out at 99.99% confidence. In the radial range 
r20oo £ r < riooo J the values of a are consistent with a 
constant and significantly larger than 1.5 (see Table 2). 
However, if only the 6 hotter ASCA clusters are fitted, 
the slope at riooo is aiooo = 1-8 ± 0.5, consistent with 3/2 
within the large errors. 

4. DISCUSSION 

4.1. The slope of the M-T relation 

Our main result is that the M-T scaling in the 1-10 
keV range is inconsistent with the self-similar prediction. 
A possible explanation for the steeper slope of the M-T 
relation, compared to the self-similar one, is preheating of 
the intracluster gas by supernova driven galactic winds be- 
fore the clusters collapse, as proposed by e.g. David et al. 
(1991), Evrard & Henry (1991), Kaiser (1991), Loewen- 
stein & Muskotzky (1996), to explain other X-ray data, 
such as the Lx — T relation and cluster elemental abun- 
dances. If supernovae release a similar amount of energy 
per unit gas mass in hot and cool clusters, the coolest 
clusters would be affected more significantly and exhibit 
a stronger shift to higher temperatures in the M-T dia- 
gram (see Figure 2) than the hotter clusters. This will 
steepen the M-T relation. The cluster formation simula- 
tions by Metzler & Evrard (1997), which include super- 
nova heating, produce a slightly steeper slope (Mcx T^-^^) 
compared to the self similar slope 1.5 in their results of 
simulations with no winds. They show that if all the wind 
energy is thermalized and retained within a virial radius, 
the temperature at masses ^ IO^^Mq may increase by 
100%, which would totally break the scaling at low tem- 
peratures. They note that in reality the effect should not 
be that dramatic since the extra energy is spent on work 
to lift the gas within the cluster potential well. 



Our results arc consistent with this explanation. Fit- 
ting only the ASCA data points (kT > 4 keV) at riooo 
with aiooo fixed to 3/2 leads to an acceptable fit with 
mooo = l-06loo9 and = 6.7 for 5 d.o.f. (see Figure 2). 
HCG62, NGC507 and NGC5044 are then 50%, 30% and 
25% hotter than what the extrapolation of the self-similar 
fit would predict for their masses. These amounts are rea- 
sonable in the supernova heating scheme, according to the 
above simulation (Metzler & Evrard 1997). 

Other work is also consistent with supernova heating 
scenario. Horner et al. (1999) examined the cluster M- 
T relation using several methods and found that spa- 
tially resolved temperature profile analysis for a sample 
of clusters with kT > 3 keV gives (T)"^/^ scaling, whereas 
their isothermal analysis sample that includes some cooler 
groups with kT > 1 keV gives a steeper slope, consis- 
tent with ours. Hjorth et al. (1998) compare tempera- 
tures and gravitational lensing masses of clusters with kT 
> 5 keV and obtain an M-T slope consistent with the 
self-similar one. The analysis of a sample of clusters with 
kT > 4 keV (Neumann & Arnaud 1999) shows that the 
self-similar slope of 3/2 that follows from their isothermal 
assumption, is consistent with their data. These obser- 
vations, together with our present results, are consistent 
with the hypothesis that at high temperatures the M-T 
scaling is self-similar, but breaks down at low tempera- 
tures (~ 1 keV). However, Horner et al. (1999) also com- 
pare temperatures with masses determined from velocity 
dispersions using the virial theorem. That dataset appears 
to be consistent with the self-similar scaling, even though 
their sample contains clusters with temperatures as low as 
~ 2 keV. The source of this disagreement is unclear (see 
4.2. for more discussion on the virial masses). 

Also supportive of energy injections, the work of Pon- 
man. Cannon & Navarro (1999) shows that cool (T < 4 
keV) clusters observed with ROSAT and GINGA have en- 
tropics higher than achievable through gravitational col- 
lapse alone which they explain by pre-heating from strong 
galactic winds. 

In the hydrostatic equilibrium scheme, since approxi- 
mately MioL oc T X /? at a given radius, for a given object 
with a certain total mass, a temperature rise due to pre- 
heating will be compensated by a shallower gas density 
profile. If heating is more prominent at small tempera- 
tures, one would then expect lower values of in cooler 
objects. This seems to be consistent with observations 
(e.g. Mohr & Evrard 1997, Vikhlinin et al. 1999, Ponman 
et al. 1999). Our sample also exhibits this behaviour. 
The average values of (3 with rms errors for the ROSAT 
and the ASCA subsample are 0.54 ± 0.09 and 0.70 ± 0.06, 
respectively. 

4.2. The normalization 

-^1000 values given by our best fit model with the ex- 
ponent as a free parameter are 2.8 and 1.4 times smaller 
than corresponding values obtained by the Evrard et al. 
(1996) n = 1 simulations at T = 1 keV and 10 keV, re- 
spectively. Note that both ASCA and ROSAT data, which 
have been obtained using independent temperature mea- 
surement methods and different instruments, give smaller 
values compared to the simulations. The above mentioned 
temperature profile analysis by Horner et al. (1999) gives 
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a normalization 40% lower than Evrard et al. (1996) at 
r'200, which is the same difference as we find between the 
normalization of our a = 3/2 fit and Evrard et al. (1996) 
£it ^1000 • The normalization of the isothermal sample re- 
sults of Neumann & Arnaud (1999) also is 30% lower than 
the Evrard et al. (1996) values at r = 0.3 r2oo- It is not 
inconceivable that the X-ray-measured masses are a factor 
of ^^2 lower than the true masses (in simulations), due, for 
example, to significant gas turbulence or magnetic fields 
(e.g. Loeb & Mao 1994) invoked to explain the differ- 
ence between the X-ray and lensing mass measurements. 
However, this seems increasingly unlikely as the cooling 
flow and temperature gradient effects appear to account 
for most of that disagreement (Allen 1998; Markevitch et 
al. 1999). 

The gravitational lensing mass analysis by Hjorth et 
al. (1998) gives a 12% lower Al — T normalization than 
Evrard et al. (1996), while Sadat, Blanchard & Oukbir 

(1998) found gravitational lensing masses in their sample 
to be 36% below the Evrard et al. (1996) scaling. Lensing 
masses are highly uncertain at present (see e.g. Hjorth et 
al. 1998). Within the errors, both above results are in 
agreement with our and other X-ray results. 

On the other hand, the virial sample in Horner et al. 

(1999) gives a normalization consistent with Evrard et al. 
(1996). There are limitations for each of the three mass 
measurement methods (virial. X-ray, lensing). For exam- 
ple, the virial masses may be inflated by inclusion of back- 
ground and foreground galaxies. The comparison is best 
done on the case-by-case basis, which is out of the scope 
of this paper. The qualitative agreement of all other re- 
sults suggests that it is the simulated values by Evrard et 
al. (1996) that may be incorrect. These simulations also 
produce too steep gas density profiles and too shallow tem- 
perature profiles compared to observations (e.g. Vikhlinin 
et al. 1999, Markevitch et al. 1998, respectively). Com- 
parison of several independent cluster formation codes in- 
dicates that gas temperature and density profiles are areas 
where different simulations disagree significantly (Frenk et 
al. 1999). For example, simulations by Bryan & Nor- 
man (1997) predict temperature profiles similar to those 
observed and should therefore produce an M-T relation 
closer to that observed. 

5. CONCLUSIONS 

We studied a sample of 9 relaxed galaxy clusters, galaxy 
groups and galaxies whose temperatures range from 1-10 
keV and for which accurate temperature profiles are avail- 
able. For the hotter subsample, the hydrostatic total mass 



profiles have been accurately determined up to r^Qo using 
gas temperature profiles measured with ASCA. For the 
cooler subsample, the mass profiles are determined from 
spatially resolved spectroscopy of the ROSAT PSPC up to 
radii ri5oo - '''looo- The mass profiles of the cool subsample 
are consistent with the "universal" NEW model, as has 
earlier been found for the hotter subsample. We derived 
the mass-temperature relation at riooo over the 1-10 keV 
temperature range as 



Minnn = 1.23 ± 0.20 10^^ h^^ M, 



\10 keVy 



1.79±0.14 



(3) 

The normalization is significantly smaller than that pre- 
dicted by simulations of Evrard et al. (1996) with = 1. 
Our relation is significantly steeper compared to the self- 
similar one (slope of 3/2). However, fitting only our ASCA 
data (kT > 4 keV) at riooo with aiooo fixed to 3/2 leads 
to an acceptable fit 



Miooo = 1-06 ± 0.09 10^5 h^o' Mq (^^^ 



1.5 



(4) 



with 2.6 fj lower normalization than in the above men- 
tioned simulations (Evrard et al. 1996). Although we can- 
not exclude a single power law slope over the whole tem- 
perature range, this behaviour is consistent with a break 
in the self-similar scaling at low temperatures, as expected 
if preheating of the intracluster gas by supernova driven 
galactic winds has taken place. Most independent M-T 
relation observations are consistent with our results and 
with the preheating scenario. The gas density slopes in 
our low temperature sample are smaller than the ones in 
the hot subsample, consistent with this supernova heating 
scheme. 

XMM and Chandra missions will be useful for extending 
the sample size. Their large effective area, excellent an- 
gular resolution and better energy resolution over a wide 
energy range will improve the accuracy of the spatially re- 
solved spectroscopy manyfold compared to the presently 
available data. 
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Table 1 
The data 



name 


''2000 


''ISOO 


''1000 


''500 


Af2000 


-^^1500 


A^iooo 


M500 


< T > 


z 




Mpc 


Mpc 


Mpc 


Mpc 


10"Mq 


10"Mq 


lO^M© 


IO^Mq 


keV 





A2029 





87 


1 


02 


1 


23 


1 


68 


4.8+1-? 


c 7+2.0 


6.61^1 


^■3-3.6 


9.10 ± 1.0 


0.0767 


A401 





92 


1 


06 


1 


27 


1 


73 










8.0 ±0.4 


0.0748 


A2256 





9f 


1 


07 


1 


28 


1 


73 


5 2+1-1 

"J-^-l.O 




7 q+1.5 

' -5-1.5 


8.8l}:« 


7.3 ±0.5 


0.058 


A3571 





9f 


1 


05 


1 


26 


1 


68 




5.61- 






6.9 ±0.2 


0.040 


A2f99 





7f 





81 





96 


1 


31 


9 9+0.5 
^■^-0.5 


9 c;+0.5 


9 0+0.6 


c:+0.7 
«5-"J-0.7 


4.8 ±0.2 


0.0299 


A496 





73 





83 





99 


1 


33 


9 c+O.S 


9 0+0.6 


q 1 +0.6 

"-••^-o.e 


q 7+0.8 
"J- ' -0.8 


4.7 ±0.2 


0.0331 


HCG62 





30 





34 





42 





57 


"■^"J-o.os 


"■^^-0.06 


'-'•^^-0.07 


n 9Q+0.09* 
'-'•^'^-0.09 


1.16 ±0.08 


0.0138 


NGC5044 





32 





37 





45 





62 


o.2ol°:°l 


r, 90 + 0.05 

'-'•^'-'-0.05 


n 97+0.05* 
^■^'-0.05 


n q£?+0.07* 
L'-'J"J_o.07 


1.09 ±0.03 


0.0087 


NGC507 





31 





35 





44 





68 




20+° °'i 

"■^^-0.04 


25+° °''^ 

^•^0-0.05 


^■^-'-0.10 


1.09 ±0.03 


0.0162 



Masses and radii are evaluated using Hq — 50 km s ^ Mpc ^. The tabulated temperatures are the observed 
ones, not redshift-corrected. Errors are given at 90% confidence. 



* Excluded from the fit 



Nevalainen J., Markevitch M and Forman W. 
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Table 2 
Fit results 





^'2000 

a = 1.5 


a free 


^1500 

a EE 1.5 


a free 


»'iooo 
a = 1.5 


a free 


a = 1.5 


a free 


n 

a 

xVd.o.f. 


33.9/8 


0.93±0.14 
1.77±0.11 
7.6/7 


0.73±0.05 
34.5/8 


1 n5+°-i^ 

1 77+0.12 
' -0.11 

7.4/7 


0.89tg:°? 
29.2/7 


-I r,q + 0.21 

1 70+0.14 

6.4/6 


1.28±0.12 
5.1/5 


1 ct;+0.50 
l.OO_0.43 
1 0^+0.51 

3.2/4 



NGC5044 HCG62 NGC507 




051015 20 25 5 1015 20 5 10 15 

Radius [arcmin] Radius [arcmin] Radius [arcmin] 



Fig. 1. — The mass profiles for the ROSAT subsample (solid line) with 90% confidence errors (dashed hue). The profiles 
are plotted in the radial range beyond the cooling flow region and up to the maximum radius of the significant surface 
brightness detection. The vertical dash-dot line shows the maximum radius of temperature measurement (in case of 
NGC507 it is equal to the maximum radius of brightness detection). The dotted line shows the total mass model with 
NFW profile for the dark matter component. 
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Fig. 2. — Temperatures (scaled to z = 0) and Afiooo values with f a errors are plotted as crosses. The best fit model to 
the ASCA data with the exponent fixed to 1.5 is shown as a dashed line. The best fit powerlaw to all data (excluding 
NGC5044, which is shown for illustration only) with exponent as a free parameter is shown as a solid line. The scaling 
law predicted by simulations (Evrard et al. 1996) is shown as a dotted line. 



